Microscopic theory of equilibrium properties of 
F/S/F trilayers with weak ferromagnets 
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The aim of this paper is to explain the non monotonic temperature dependence of the self- 
consistent superconducting gap of ferromagnet/superconductor/ferromagnet (F/S/F) trilayers with 
weak ferromagnets in the parallel alignment (equivalent to F/S bilayers). We show that this is 
due to Andreev bound states that compete with the formation of a minigap. Using a recursive 
algorithm we discuss in detail the roles of various parameters (thicknesses of the superconductor and 
ferromagnets, relative spin orientation of the ferromagnets, exchange field, temperature, disorder, 
interface transparencies) . 
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I. INTRODUCTION 

In conventional superconductivity the attractive inter- 
action mediated by phonons binds electrons into Cooper 
pairs that condense in the BCS ground state with a zero 
temperature gap A to the first quasiparticle excitations 1 . 
In ferromagnetism electron interactions generate a spin 
symmetry breaking that can be described by the Stoner 
model in which electrons subject to an exchange field h ex 
acquire a Zeeman energy. 

Many physical phenomena are involved at the inter- 
faces between superconductors (Ss) and ferromagnets 
(Fs). For instance it was shown in the early 1970's that 
the Fermi surface spin polarization of a ferromagnetic 
metal could be measured by spin-resolved tunneling be- 
tween a ferromagnet and a superconducting film in the 
presence of Zeeman splitting^.. The Fermi surface spin 
polarization was measured more recently^ by Andreev 
reflection at F/S interfaces^ with highly transparent in- 
terfaces, not in thin film geometries. The non equilib- 
rium spin population in the ferromagnet plays also a role 
in Andreev reflection at F/S interfaces*^* 8 -. Andreev re- 
flection with Zeeman splitting in a thin film geometry 
could be another way to probe the Fermi surface spin 
polarization 9 . 

These experiments can be interpreted without dis- 
cussing the proximity effect (the pair correlations in- 
duced in the normal metal or ferromagnet) that has fo- 
cused an important interest recently. It is well estab- 
lished both theoretically and experimentally that the 
pair amplitude induced in a ferromagnetic metal oscil- 
lates in space and can become negativei2*iiii2*iilii4 ; giv- 
ing rise to the 7r-coupling. The 7r-coupling generates 
oscillations of the critical temperature of F/S multi- 
layers as a function of the thickness of the ferromag- 
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netic layersi£*i§*i£*i2*i**V The bad quality of interfaces 
may however also play a role in these experiments 2 ^* 2 ^. 
Direct evidences of the 7r-coupling have been obtained 
recently 22 i 23 i 24 ' 25 . The proximity effect at F/S interfaces 
is related to the formation of the so-called Andreev bound 
states. Andreev bound states were first discussed by 
de Gennes and Saint-James 2 ^ and Andree\*2£ for normal 
metal (N)/S interfaces. Several theoretical investigations 
of Andreev bound states at F/S interfaces have been pre- 
sented recently 21 - 28-29 as well as numerical investigations 
of the proximity effect at F/S interfaces, based on simu- 
lations of the Bogoliubov-de Gennes equations^. 

The interest in the inverse proximity effect (the prop- 
erties of the superconductor) at F/S interfaces dates back 
to the 1960's3ii22i2ii where it was shown that with insulat- 
ing ferromagnets an exchange field is induced in the su- 
perconducting electrode of a F/S/F trilayer in the paral- 
lel alignment. The theoretical prediction by de Gennes^i 
was further confirmed by experiments in the late 1960's 
with metallic 32 and insulating 33 ferromagnets. Recent 
experiments were carried out with metallic ferromagnets, 
which confirmed the effect on the critical temperature 3 ^. 
It was shown recently 3 ^*! that an exchange field in the 
superconductor exists also with metallic ferromagnets 
but the sign of the magnetization in the superconduc- 
tor is opposite to the magnetization in the ferromagnet. 
An exchange field in a superconductor is pair breaking. 
As a consequence the superconducting transition tem- 
perature in the parallel alignment is smaller than in the 
antiparallel alignment. 

With metallic ferromagnets the proximity effect may 
influence the value of the superconducting gap and tran- 
sition temperature. It was realized recently^* 3 ^*^^* 3 ^ 4 ^ 
that the proximity effect in F/S/F trilayers is quite 
special since under some conditions the zero temper- 
ature superconducting gap in the parallel alignment 
can be larger than in the antiparallel alignment. This 
was established 3 ^*^ within a model of multitermi- 
nal hybrid structure originally proposed for transport 
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and is related to spatially sepa- 



rated superconducting correlations among the two ferro- 
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magnets. The same behavior was found within a model 
of F/S/F trilayer with atomic thickness and half- metal 
ferromagnets^, which was finally extended to Stoner 
ferromagnets^. 



Other physical effects take place with weak ferromag- 
nets for which the exchange field is smaller than the su- 
perconducting gap, as indicated by the reentrance of the 
critical temperature of F/S superlattices as a function 
of the exchange field^S^i, and by the reentrance of the 
critical temperature of F/S bilayers as a function of the 
thickness of the ferromagneteSi224ii. The full tempera- 
ture dependence of the self-consistent superconducting 
gap was calculated in Ref. and it was shown that for 
F/S bilayers or F/S/F trilayers with atomic thickness the 
superconducting gap can have a non monotonic tempera- 
ture dependence: within a given range of parameters the 
superconducting gap first increases as temperature T is 
reduced, reaches a maximum, and decreases to zero as T 
is further reduced. Within a narrow range of interface 
transparencies there exists also a reentrant behavior of 
the superconducting gap at low temperature^. 
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FIG. 1: (a) Schematic representation of the F/S/F trilayer. 
D and H are sent to to infinity (infinite planar geometry). 
The ferromagnets (superconductor) are made of Lf (Ls) lay- 
ers stacked in the lateral (z) direction, (b) Schematic repre- 
sentation of the 2D F/S interface model-fii considered in sec- 
tion ^^^| d! is sent to infinity. 



The purpose of this article is to show that this behav- 
ior is related to the formation of Andreev bound states 
that can compete with the formation of a superconduct- 
ing minigap. We find a systematic correlation between 
Andreev bound states at the Fermi level and a reduction 
of the low temperature self-consistent superconducting 
gap, which constitutes the main result of this article. 



The article is organized as follows. Preliminaries are 
given in section ^ in which we discuss the models and 
the Green's functions formalism. A recursive algorithm, 
more efficient that the direct inversion of the Dyson 
matris*^, is presented in section IIIII Different regimes 
can be obtained depending on how the lateral dimensions 
Ls of the superconductor and Lp of the ferromagnets 
compare respectively (i) to the zero temperature lateral 
superconducting coherence length £g — 2i ao/7rA2&, 
where to is the lateral hoping amplitude, and ao is equal 
to the interatomic distance in the tranverse direction; 
and (ii) to the zero temperature lateral exchange length 
£p = 2toao/7r/i e ^2&. We present a detailed investiga- 
tion of the different regimes 

L F ,L 5 <Cs L) f section !!^ : 
Ls < 4 ±) and L F < ^ (section Q; L s < ^ and 
Lp > (section fVTl . We cannot carry out a system- 
atic study of the regime Ls ~ , Lp > £p which is 
too demanding from a computational point of view. In 
section lYl!! we provide a comparison with other models 
proposed recently to describe Andreev bound states in 
F/S IrybridsSSiSi. Concluding remarks are given in sec- 
tion E3H1 
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FIG. 2: Cut in the (z,y) plane of the F/S/F trilayer with 
(Ls/ ao, I/f/oo) = (3,4), where ao is the interplane spacing. 
The ferromagnets and the superconductor are infinite in the 
x and y directions. The left ferromagnet ends at the layers 
a' and a. The right ferromagnet ends at the layers b and b' . 
The superconductor ends at layers a and /3. 



II. PRELIMINARIES 

A. The models 

We suppose that the two ferromagnets of the F/S/F 
trilayer have the same thickness Lp. We note Ls the 
thickness of the superconductor (see Fig. The three 
electrodes are made of 2D planes stacked along the z axis 
(see Fig- EJ - Each superconducting plane is described by 
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the BCS Hamiltonian 



B. Green's functions 



(1) 



1. Zero temperature Green's functions 



+ A 12 ( c tT c -k,i + c k,l c -k,T) , 



where the wave vector k is parallel to the 2D layer. The 
free electron dispersion relation within one layer is c(k) = 
H 2 k 2 /2m. The ferromagnets are described by the Stoner 
model 



H 



Stoner — 1^. e (k) c k.g C k,° 



(2) 



k,cr 



T Ck,T - c ka c ka 



where h cx is the exchange field. The coupling between 
two planes belonging to the same superconducting or fer- 
romagnetic electrode is given by 



W F _ F = W S -S = to J2 



(3) 

whereas the coupling between the ferromagnets and the 
superconductor is given by 



(4) 



where the summation runs over all sites at the inter- 
face (xi corresponds to a site on the left side and Xp 
is the corresponding site on the right side). The pa- 
rameters are such that the coherence length within one 
layer £|P = Kvp /A and the tranverse coherence length 
= 2toaQ/irA are larger than the width L$ of the su- 
perconductor. The model should stricktly speaking apply 
to ballistic systems whereas real samples are usually in 
the diffusive regime. However we can include disorder 
in one particular case (the F/S/F trilayer with atomic 
thickness, see section HV C|) . In this case the qualitative 
behavior is robust against increasing disorder. Fabry- 
Perot resonances generate parity effects for small values 
of Ls and Lp in the ballistic model. These parity effects 
are not expected to occur in the diffusive regime, but 
do not occur either in our simulations if Ls and Lp are 
sufficiently large. By increasing Lp and Lg we obtain 
a cross-over between the regime Lp,L$ < •Vf an d the 
regime Xp < Lp,Ls <C £s • We find the same quali- 
tative effects in the two regimes. Moreover we obtain in 
section^a non monotonic variation of the self-consistent 
superconducting gap as a function of the exchange field. 
This is compatible with the non monotonic variation of 
the critical temperature as a function of the exchange 
g e i ( j38 j 48 J 49 obtained in the context of linearized Usadel 
equations for disordered conductors. The compatibility 
between the two behaviors indicates the validity of our 
approach. 



The Green's functions of an isolated superconduc- 
tor can be gathered in a 4 x 4 matrix in the spin ® 
Nambu representation but in the absence of non collincar 
magnetizationsMi£2*£iiii£ the quantization axis can be 
chosen parallel to the exchange field so that the 4x4 
Green's functions reduce to two separate 2x2 matrices, 
one in each spin sector. For practical purpose we work 
in the spin-up sector. The Green's function is given by 

ffx,y(t,t') = -i (5) 

(T t fcx, T (t),c+ T (t')l) (T t (cx, T (t),c y a (*'))> 

( T *(<i(0:<t(t')J) (^(c+^C^ft'))) 

where x and y are two arbitrary sites in the supercon- 
ductor and T t is the usual T-productS. The "11" com- 
ponent describes the propagation of a spin-up electron, 
the "22" component describes the propagation of a spin- 
down hole and the "12" and "21" components describe 
superconducting correlations. After Fourier transforming 
we obtain a standard expression for the different elements 
of the Green's function 53 : 



J a, a VS : 



E k +irj oj - 
A 



E k 



ii] 



uj = - 



[uj - E k + it]} [uj + E k - irj\ ' 



(6) 
(7) 



The variable £ is related to the kinetic energy: £, k = 
h 2 k 2 / (2rn)—ep where ep = h 2 k 2 F /2m is the Fermi energy. 
E k = ^/A 2 + ^ is the quasiparticle energy and u\ = (1 + 
£, k /E)/2 and v 2 k = (1 - £k/E)/2 are the BCS coherence 
factors. 

The Green's function of a spin-up ferromagnet is di- 
agonal in Nambu space. The "11" component is given 
by 

W) = \u> - I + h cx + ii] sgn(£ - h cx )} ' (8) 

and the Green's functions of a spin-down ferromagnet are 
obtained by changing h ox into —h cx . 

The Green's functions G x . y of the connected trilayer 
are given by the Dyson equation G = g+g®T,®G, where 
in a compact notation S is the self-energy corresponding 
to the tunnel Hamiltonian (0J and ® corresponds to a 
summation over spatial variables and a convolution over 
time variables. We look for non perturbative solutions 
of the Dyson equations suitable for describing Andreev 
bound states. 



2. Finite temperature Green's functions 

Finite temperature Green's functions are obtained 
through the analytic continuation ui — > iuj and by sum- 
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ming over the Matsubara frequencies u n = (2n + 1)tvT 
where T is the temperature^. The superconducting gap 
is determined from the BCS self-consistency equation^ 



d 2 k 



(9) 



where A is the strength of the attractive electron-electron 
interaction. To evaluate we change variable to £ = 
H 2 k 2 /(2m) — ep and restrict the integral to |£| < u>n- 
To avoid introducing new parameters we use ujd = £f = 
h 2 k F /2m. We note Ao the superconducting gap of an 
isolated 2D layer. All energy scales will be compared to 



III. RECURSIVE ALGORITHM FOR THE 
F/S/F TRILAYER 

A. Green's functions of an isolated electrode 

We aim to describe F/S/F trilayers with a fi- 
nite thickness of all electrodes. In this respect the 
(Ls I 'ao, Lf I 'ao) = (1,1) trilayer is just viewed as a toy- 
model: we use a mean field approach that does not in- 
corporate the phase fluctuations of the order parameter 
and we do not consider possible instabilities such as spin 
or charge density wave. 

The ferromagnetic and superconducting electrodes of 
the F/S/F trilayer consist of a finite number of layers 
stacked along the z axis and labeled from 1 to I (see 
Fig. El- Two consecutive layers n and n+1 are coupled by 
a tunnel amplitude t n . We use Green's functions that are 
parametrized by the wave vector in the (x, y) direction 
and by the spatial coordinate in the lateral direction. 

We note h[ L j the Green's functions of the system of L 
layers in a given electrode and gj,, the Green's function 

of the isolated layer number i. The Green's function h^ L 
can be calculated recursively through a matrix continued 
fraction: 



L L.L 



f - S iiL- 1 ) $ 

1 — 9L,LtL-l'l L _ 1 L _ 1 lL-l 



9L,L- 



(10) 



The local Green's functions of the system of L stacked 
layers are obtained through 



where and h^l are calculated recursively through 
the relations — h^\tL-\h^Zi\ and h\ L ^ — 

h^ L L ^tL-ih^ L . The computation time required to ob- 
tain h^\ is proportional to L whereas it is proportional 
to L 2 if one wants to calculate all the h[ L ) . 



B. Green's functions of the connected trilayer 

We note G the Green's functions of the connected 
F/S/F trilayer (see Fig. |2J). Layers a and a are con- 
nected by a tunnel amplitude t a , a = t a>a and layers b 
and (3 are connected by a tunnel amplitude tb t p = tpj- 
We use the notation t = t a ^ a = tb,p and denote by to the 
tunnel amplitude within each layer (see Eq.|3|and Eq.0J. 

The definition of the tunnel Hamiltonian that we use 
for the F/S/F trilayer is slightly different from the con- 
ventional definition given by ijljl: 

WfSF = —7= ( c x,y,a,a- c x,y,a,<j + c x ,y ,a ,a c x ,y ,a ,a\^2) 

* x.er 

+ ~~7k {^x,y,p,a C x,y,b,rr + c t,y,b,<7 C x,y,P,<yj ■ 

* x.cr 

With the factors \[2 the F/S/F trilayer with Ls/ao = 1 
in the parallel alignment is equivalent to the F/S bilayer 
with Ls/ao — 1 and a tunnel amplitude equal to t. The 
correspondence between the bi and trilayer is useful for 
checking the numerical simulations. For Ls/a > 2 we 
carry out the simulations of the F/S/F trilayer but we 
note that in the parallel alignment the qualitative predic- 
tions are valid also for F/S bilayers as long as the thick- 
ness of the superconductor is smaller than the coherence 
length. 

We note K a , a = h a , a i a , a h a . a i a , a , Kp t p = 

hp,ptp,bhb,btb,p, Ka,j3 = h a ,ptpfihb,btb,Pi Kp,u = 

h/3,ata,ah a ,ata,a- The Green's function G a>a is given by 



G 



I-K n 



K, 



a,/3 









Kf3,a 



(13) 



X h a ^ a + K a fi I 

and Gp iCC is given by 







i - iW 


hp, a 



G 



p,a 



I-K, 



0,0 



-1 r 



P,a 



Kf3,aG 



(14) 



The Green's function G a ,p is deduced from Gp ta through 
the relation G^'^ 2 = G^ 2 ^" 1 , where t\ and are the 
Nambu indexes. Gp t p is given by 

l 



G 



0,0 



K 



0,0 



hp, [3 + Kp, a G a ,p 



(15) 



We deduce the values of G a , a , G&,6, G a ,b and Gb, a as well 
as Gi^ in the superconductor: 

G%^i hi i -f- hi j0 ,t aja G a ^ a t a ^ a h a ^i -\- hi^ a t a a G a ,btb,php 

+ hi : ptp v bGb, a t a ,aha,i + hi t ptp t bGb,bib,php,i. (16) 

To obtain the pair amplitude in the ferromagnets and 
superconductor we first calculate recursively the h's and 
next evaluate G\^- The evaluation of the self-consistent 
superconducting gap is done either by dichotomy if 
Ls/ao = 1, 2 or by iterations of the self-consistency equa- 
tion © if Ls/ao > 3. 
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FIG. 3: Temperature dependence of the self-consistent super- 
conducting gap for the F/S/F trilayer in the parallel align- 
ment with {L s /a ,L F /a ) = (1,1) (a), (1,2) (b), (1,3) (c), 
(1,4) (d), with weak ferromagnets (/i GX /Ao = 0.83). Lf 
and Ls are small compared to £g = 2toao/7rAo = 17.7ao, 

= 2toao/irh cx = 21.2ao and Af = 6.28ao. We use 
Aq/ef = 0.014 and t /e F = 0.4. 




IV. 



F/S/F TRILAYERS WITH L S ,L F < I^.Cf" 3 



FIG. 4: Energy dependence of the spin-up LDOS (in arbitrary 
units) in the superconducting layer of the F/S/F trilayer in 
the parallel alignment with (Ls/ao, Lp/ao) = (1,3). The 
superconducting gap, equal to Ao, is not self-consistent, (a) 
and (b) correspond to h ex /Ao = , £p = oo (N/S interface). 

(c) and (d) correspond to /i cx /A = 0.83, = 21.2a (F/S 
interface with weak ferromagnets). (e) and (f) correspond 

to h cx /Ao = 13.9, — 1.3ao (F/S interface with strong 

ferromagnets). We use Ao/ef = 0.014. to/ef = 0.4, rj/Ao = 
8.3 x 10" 3 , A F = 6.28a . 



In this section we consider the regime Ls,Lp <C 
Cs^, i establish a connection between the LDOS and 
the self-consistent superconducting gap, and show that 
the regime h ex / Aq ~ 1 is characterized by Andreev 
bound states competing with the formation of a mini- 
gap. 



nant case (see Figs. 13(a) and (c)). This is because the 
effective coupling between the superconductor and ferro- 
magnets is much weaker in the off-resonant case due to 
the differences in the density of states. 



A. Self-consistent superconducting gap for weak 
ferromagnets in the parallel alignment 

1. Role of the thicknesses of the electrodes 

We have shown on Fig. the temperature dependence 
of the superconducting gap for values of (Ls,Lp) such 

that Ls, Lp < i { s ^\^p ] . We used L s /a = 1 on Fig. El 
but similar results were obtained with Ls/ao — 2. As 
expected from Appendix ^ the breakdown of supercon- 
ductivity in the resonant F/S/F trilayer in the parallel 
alignment resembles the case (Ls/ao, Lp/ao) — (1, l) 40 : 
as temperature is reduced the superconducting gap in- 
creases, reaches a maximum and decreases to zero. We 
obtain a reentrant behavior in a narrow range of inter- 
face transparencies (not shown on Fig. [SJ. By contrast 
we obtain a monotonic behavior for off-resonant values 
of (Ls/ao, Lp/ao), a s expected from Appendix lAl 

The values of the tunnel amplitude needed to destroy 
superconductivity for off- resonant trilayers (see Figs. |3J- 
(b) and (d)) is almost ten times larger than in the reso- 



2. Role of the exchange field 

We repeated the simulations for the F/S/F trilayer 
with (Ls/ao, Lp/ao) — (2,2) but with different values 
of h cx /A . For the smallest value of h cx / A (h cx /A — 
0.28) we obtained a monotonic decrease of the supercon- 
ducting gap as a function of temperature for all values of 
the tunnel amplitude. For h ex /A = 0.56, 0.83, 1.11, 1.39 
we obtained a non monotonic temperature dependence of 
the superconducting gap similar to Fig.EJ(a) and Fig.|3 
(b) . The non monotonic variation of the superconducting 
gap occurs typically for h cx being a fraction of Ao up to 
values slightly above A . 



B. Relation between the local density of states and 
the self-consistent superconducting gap 

1. Local density of states in the parallel alignment 

The spin- up LDOS p-\(u>) in the superconductor is 
shown on Fig. 0] for the F/S/F trilayer in the paral- 
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lei alignment, for (Ls/a ,Lp/a ) = (1,3). We ob- 
tained similar results for (L s /a , L F /a ) — (2,2). The 
spin-down LDOS pi(uj) is obtained through the relation 
P[{oj) = pi(—u>). The zero temperature superconduct- 
ing gap A is fixed to the BCS value A of an isolated 
superconducting layer. 

The case of a N/S interface is shown on Fig. |3J(a) and 
Fig. lU-(b). The LDOS is symmetric with respect to a 
change of sign in energy, as expected in the absence of 
an exchange field. There exists a minigap at the Fermi 
energy so that superconductivity is robust against in- 
creasing the tunnel amplitude t. Increasing t gives rise 
to pairs of Andreev bound states at opposite energies (see 
Fig.|IJ-(a)). Each peak corresponds to a miniband, which 
is visible on Fig. 0}(b) obtained with larger values of t. 
The formation of the minibands is related to the infinite 
planar geometry: the layers are infinite in the x and y 
directions (see Fig. ^ and there is thus a degeneracy as- 
sociated to the position of the Andreev bound state in 
the (x, y) plane. The bound state can delocalize in the 
(x, y) plane and thus acquire a dispersion which would 
not occur if the dimensions H and D (see Fig. ^| were 
small compared to the BCS coherence length, a situation 
considered in Ref. |2^. 

For weak ferromagnets with t, h cx and A having the 
same order of magnitude (panel (c) on Fig. |IJl we obtain 
one Andreev bound state miniband inside the gap and 
one resonant scattering state above the gap. The An- 
dreev bound state miniband moves to the Fermi energy 
as the tunnel amplitude t increases (see Fig. IU-(c)) and 
appears at a positive energy for larger values of t (see 
Fig.H(d)). 

For larger values of the exchange field the Andreev 
bound state miniband disappears from the LDOS. For 
/iex/Ao 3> 1 the induced exchange field is opposite to 
the magnetization in the ferromagnets*^*^. 

There are thus two qualitatively different depairing 
mechanisms for strong and weak ferromagnets. For 
strong ferromagnets with h ex / Ao S> 1 (see Figs. |U(e) 
and |U-(f)) the breakdown of superconductivity is due 
mainly to Zeeman splitting. The case of very weak ferro- 
magnets (h cx /A <C 1) resembles that case h ex /A = 0: 
non perturbative Andreev bound states are generated in 
the superconducting gap approximately at opposite en- 
ergies. In the case h cx /Ao ~ 1 the Andreev bound state 
miniband can be at zero energy, therefore destructing the 
minigap. 



2. Finite temperature local density of states 

The spin-up LDOS at a finite temperature T is re- 
lated the conductance of a scanning tunneling microscope 
(STM) in which the tip is made of a half-metal ferromag- 
net. The finite temperature current at an arbitrary volt- 
age is obtained through Keldysh formalism. The finite 
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FIG. 5: Energy dependence of the finite temperature 
spin-up LDOS (in arbitrary units) in the superconducting 
layer of the F/S/F trilayer in the parallel alignment with 
(Ls /ao, Lf /ao) = (1,3). The zero temperature LDOS cor- 
responds to Fig.El(c) witn t/Ao = 1.4. 



temperature LDOS is found to be 

p(u>) 



Pt(u) — / du>' 



4Tcosh 2 (^)' 



(17) 



where p(oj) is equal to the zero temperature LDOS. 

We calculated the finite temperature LDOS of the 
(I/s/ao, Lp/ag) = (1, 3) trilayer in the parallel alignment 
(see Fig. |SJ). Increasing temperature tends to reduce the 
intensity of the Andreev bound state at the Fermi en- 
ergy on Fig.QJ-(c) so that the LDOS at the Fermi energy 
decreases if T increases. For relatively large values of 
T the peak structure has almost disappeared from the 
LDOS but there remains a minimum associated to the 
superconducting gap that is significantly filled. As T in- 
creases there is thus a first cross-over where the peak at 
the Fermi energy disappears and a superconducting mini- 
gap is restored, and a second cross-over where the super- 
conducting gap disappears. The first cross over occurs 
at a temperature equal to the bandwidth of the Andreev 
bound state miniband and the second cross-over occurs 
at a temperature comparable to the zero temperature 
superconducting gap. The finite temperature LDOS is 
thus in a qualitative agreement with the non monotonic 
self-consistent superconducting gap. 

Reentrance obtained in a narrow range of interface 
transparencies^ can also be explained by this qualitative 
picture: if the Andreev bound state miniband is narrow 
and located at a slightly positive energy like on Fig. 0J- 
(c) then the LDOS at the Fermi energy is vanishingly 
small at T = since the Fermi energy is not in the An- 
dreev bound state miniband. By increasing T the width 
of the Andreev bound state miniband increases so that 
the density of states at the Fermi energy increases. By 
further increasing T the intensity of the Andreev bound 
state miniband is reduced, and the density of states at the 
Fermi energy decreases. This behavior is compatible with 
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FIG. 6: Energy dependence of the spin-up LDOS (in arbi- 
trary units) in the superconducting layer of the F/S/F trilayer 
in the antiparallel alignment with (Ls/ao, Lf/clq) = (1,3). 
The superconducting gap, equal to An, is not self-consistent, 
(a) and (b) correspond to ft cx /An = 0.83, £p = 21.1oo 
(weak ferromagnets) . (c) corresponds to /i cx /An = 2.78, 
= 6.4an (strong ferromagnets). We used Ao/ef = 0.014, 
t /e F = 0.4, r)/A = 8.3 x 10" 3 and X F = 6.28o . 

a reentrant behavior of the self-consistent superconduct- 
ing gap. The correlation between the low temperature 
superconducting gap and the LDOS at the Fermi energy 
is further established in section IVT1 



3. Local density of states in the parallel alignment 

The spin-up LDOS in the antiparallel alignment is 
shown on Fig. H3 The LDOS is symmetric with respect 
to a change of sign of energy, but not equivalent to the 
LDOS of a N/S interface (see Fig. [1(a) and Fig.H(b)). 
There exists a well-defined minimum at the Fermi energy 
corresponding to the superconducting minigap. The en- 
ergy dependence of the LDOS in the antiparallel align- 
ment shows that superconductivity is stronger than in 
the parallel alignment, both for weak and strong ferro- 
magnets. This is expected because of the exchange field 
induced in the superconductor in the antiparallel align- 
ment is reduced compared to the parallel alignments^!. 
The self-consistent superconducting gap has a monotonic 
temperature dependence in this case which is because a 
well-defined minigap is obtained in the LDOS. 



C. Role of disorder for weak ferromagnets and for 
the F/S/F trilayer with atomic thickness 

Disorder plays an important role in the proximity ef- 
fect. Diffusive N/S interfaces are described by quasi- 
classical theory^. A small disorder can be incorporated 
in our description based on microscopic Green's func- 
tions like in Ref. (see Appendix . The strength 
of disorder in the superconductor is characterized by 
Sg = yJuaU^/ Ao, where n a is the concentration of im- 
purities and u a is the scattering potential, and we use 
a similar parameter Sp^ to characterize disorder in the 




0.25 0.5 0.75 0.25 0.5 0.75 1 



T/A 

FIG. 7: Temperature dependence of the self-consistent su- 
perconducting gap of the F/S/F trilayer in the parallel align- 
ment with (Ls/ao, Lp/ao) = (1, 1) in the presence of disorder 
in the superconducting layer. Without vertex corrections we 
use h eK /A = 0.83, 4°' = 6.2 (a), 4°' = 7.6 (b), 4°' = 9-8 
(c). (d) corresponds to 4°' — 9-8 with vertex corrections. In 
all cases we use 8p = 0. 

ferromagnets. We obtain a significant effect of disorder 
for relatively large values of and 5jP . 

1. Disorder in the superconducting and ferromagnetic 
layers 

The temperature dependence of the self-consistent su- 
perconducting gap in the presence of disorder in the su- 
perconductor is shown on Fig. [7| for (Ls/a,Q, Lp/ao) = 
(1,1) in the parallel alignment. In the absence of disor- 
der Fig.[7|corresponds to Fig. 131(a). The effect of disorder 
is to reduce the effect of the tunnel amplitude^ so that 
if disorder increases a larger value of t is needed to de- 
stroy superconductivity. Comparing Fig. [21 and Fig. we 
see that the variations of A(T) are affected by a weak 
disorder especially if the F/S/F trilayer in the parallel 
alignment is close to the breakdown of superconductiv- 
ity. In this case the dimensionless parameter controlling 
the strength of disorder is 5s = S^Aq/A(T) which can 
be much larger than 5^ . 



2. Vertex corrections 

There exist two perturbative series: one in the hop- 
ping amplitude t and the other in the disorder scattering 
potential u. Vertex corrections arise from diagrams that 
mix the two series (see Appendix IB 21 . The tempera- 
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FIG. 8: (a) Variation of the self-consistent superconduct- 
ing gap calculated at T/Ao = 0.14 in the parallel (P) align- 
ment (open symbols) and in the antiparallel (AP) align- 
ment (filled symbols) as a function of h ax /Ao- (b) Vari- 
ation of the self-consistent superconducting gap calculated 
with hsx/Ao = 0.56 in the parallel alignment, as a function of 
T/Aq. We used in both cases {Lp/ao, L s /a ) = (11, 11). The 
Fermi wave length is Xf = 6.28ao. We use t/Ao = 5.6 (□ in 
the P alignment and ■ in the AP alignment), t/Ao = 8.3 (o in 
the P alignment and • in the AP alignment), and t/Ao = 11.1 
(A in the P alignment and A in the AP alignment). 



ture dependence of the self-consistent superconducting 
gap with the vertex corrections is shown on Fig. 0-(d). 
The critical temperature is larger if vertex corrections 
are included. We obtain a non monotonic temperature 
dependence of the superconducting gap even in the pres- 
ence of vertex corrections. The role of vertex corrections 
increases if t increases since the vertex correction term is 
proportional to n a t it„. 



V. FINITE THICKNESS IN THE 
SUPERCONDUCTOR (L s < ^ AND L F < ^ F r> ) 

In this section we include a finite thickness in the su- 
perconductor for Lp <C £p '. For weak ferromagnets in 
the regime h ex / Aq ~ 1 we obtain non monotonic temper- 
ature dependences of the self-consistent superconducting 
gap in the parallel alignment, therefore confirming sec- 
tion |IV| 

We have shown on Fig. 0-(a) the variation of the self- 
consistent superconducting gap in the middle of the su- 
perconductor as a function of the reduced exchange field 
/lex/Ao- As expected from the variation of the critical 
temperature 3 ^*^^ we obtain a minimum if the exchange 
field is comparable to the superconducting gap. Moreover 
we obtain a minimum also in the antiparallel alignment 
that is related to the minigap formed in between Andreev 
bound states at opposite energies (see section IIV B 31) . 
The full temperature dependence of the self-consistent 
superconducting gap in the parallel alignment is shown 

on Fig.|SJ.(b) for Ls/ao = 11 (comparable to £g V a o) an d 
Lp/ao = 11 (much smaller than /da). We obtain a 
maximum in the variation of A(T)/Aq in the parallel 
alignment. We carried out the same simulation in the 
antiparallel alignment and found a monotonic temper- 



ature dependence of the self-consistent superconducting 
gap (not shown on Fig- - 



VI. FINITE THICKNESS IN THE 
FERROMAGNETS (L s < ^ AND L F > ^) 

We discuss now the regime 

Ls<& ) and Lp>^ ] . 
We find Andreev bound states at the Fermi energy in 
the parallel alignment for strong ferromagnets, correlated 
with non monotonic temperature dependences of the self- 
consistent superconducting gap. 

The regime Lp > £p is characterized by oscillations 
of the self-consistent superconducting gap, critical tem- 
perature and pair amplitude as a function of Lp. We 
obtain bound states within the superconducting gap and 
resonant scattering states outside the superconducting 
gap. We calculated systematically the LDOS and the 
self-consistent superconducting gap for Ls/ao = 1,2 and 
Lp/ ao = 1, 100, in the parallel and antiparallel align- 
ments (see Fig. |5J). The systematic calculation of the 
finite temperature LDOS given by Eq. (|17|) is too de- 
manding from a computational point of view. Instead 
we calculate the LDOS at zero temperature with r\ = T. 
The consistency between the two calculations was verified 
in a few cases. Depending on the interface transparen- 
cies A(Lp/ao)/ 'Ao tends either to a finite value or to zero 
in the limit Lp/ao — > +oo. We concentrate on the first 
case only. For 1 < Lp/ao < 40 we see on Fig. [3] that 
on average the superconducting gap is smaller when the 
LDOS at the Fermi energy is larger. For Ls/ao = 1 and 
i/A = 2.8 (see Fig. H(b)) we obtain A(T)/A = 
for four values of Lp/ao (Lp/ao — 4,6,11,13). For 
three other values of Lp/ao (Lp/ao — 16, 23, 25) we ob- 
tain a non monotonic variation of A(T)/Ao. For two 
other values of Lp/ao (Lp/ao — 32 and Lp/ao — 35) 
A(T)/Ao is monotonic but far from the BCS variation. 
These nine values of Lp/ao with an anomalous temper- 
ature dependence of the self-consistent superconducting 
gap have large values of the LDOS at the Fermi energy 
(p(u = 0) > 0.045 on Fig. 13(b), in arbitrary units). Like 
in section IIVB II Andreev bound states near the Fermi 
energy correlate with unconventional temperature depen- 
dences of the superconducting gap. 

For 40 < Lp/ao < 100 corresponding to Lf > £p , 
we obtain a cloud of points with a small dispersion and 
with no correlation between the superconducting gap and 
the LDOS at the Fermi energy. This corresponds to the 
cross-over to the F/S/F trilayer with bulk ferromagnets. 

We obtain bound states within the superconducting 
gap and unconventional temperature dependences of the 
superconducting gap with strong ferromagnets in the par- 
allel alignment. By contrast with strong ferromagnets we 
obtain a conventional temperature dependence of the su- 
perconducting gap for the F/S/F trilayer with smaller 
values of (Ls/ao, Lp/ao)^- 

We carried out the same simulation in the antiparallel 
alignment and found that A(T)/Aq is close to the BCS 
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FIG. 9: Correlation between A/Ao and the LDOS at the 
Fermi energy (in arbitrary units) for Ls/a0 — 1 and Lp/ao 
between 1 and 40. The temperature is T/Ao = 0.14. In the 
calculation of the LDOS the superconducting gap is not self- 
consistent. We use strong ferromagnets with /i cx /Ao = 13.9 
and £p — 1.3ao- We use Ls/ao = 1 and t/Ao = 2.1 (a); 
Ls/ao = 1 and t/Ao = 2.8 (b); L s /a = 2 and t/Ao = 3.5 
(c); Ls/ao = 2 and t/A = 4.2 (d); Ls/ao = 1 and t/Ao = 
2.8 (e); L s /ao = 1 and t/A = 3.3 (f). (a), (b), (c), (d) 
correspond to the parallel alignment, (e) and (f) correspond 
to the antiparallel alignment. The solid line in (b) corresponds 
to p{u = 0) = 0.045 (see text). 



temperature dependence of the superconducting gap for 
all values of Ls/ao between 1 and 100. A larger density 
of states at the Fermi energy correlates with smaller val- 
ues of the self-consistent superconducting gap like in the 
parallel alignment (see Fig. EJ(e) and (f)). Because of 
the minigap that can be formed since the bound states 
are at opposite energies (see Fig. |H} we do not obtain 
the points with large values of the density of states at 
the Fermi energy like in the parallel alignment. The ab- 
sence of large values of the LDOS at the Fermi energy 
is related to the fact that the Andreev levels in the an- 
tiparallel alignment do not cross the Fermi energy as the 
interface transparency is increased. 



VII. COMPARISON WITH OTHER MODELS 



The discussion in the preceding sections was restricted 
to the infinite planar geometry. Now we discuss Andreev 
bound states in other geometries without imposing self- 
consistency on the superconducting gap. In section lYlI Al 
we improve the discussion of a model proposed by re- 
cently Vecino et alm^ in which a ferromagnetic chain is 
connected to a superconductor—. In section IVII Bl we 
compare the LDOS to a model discussed recently by 
Cserti et al^. 
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FIG. 10: Reduced energy E/ A of the Andreev bound states as 
a function of the number of sites Lp/ao in the ferromagnetic 
chain for weak ferromagnets Qi ex /A = 0.5). We used the 
parameters A/t = A /t F = 0.01, t/t s = 0.01 , L s = 10 5 a . 
The entire spectrum is kept in (a). Only the levels with a 
residue larger than 2 are kept in (b). The superconducting 
gap is not self-consistent. 



A one dimensional model 



We suppose that the superconductor and ferromagnets 
are described by ID chains with open boundary condi- 
tions with Ls/ao sites in the superconductor and Lp/ao 
sites in the ferromagnet. We denote by ts (£f) the hop- 
ping amplitude in the superconductor (ferromagnet) and 
we use ts = tp. The energy level spacing in the super- 
conductor is much smaller than the superconducting gap 
A. 

We have shown on Fig. ^| the evolution of the energy 
of Andreev bound states as a function of the length of 
the ferromagnetic chain obtained with the algorithm pre- 
sented in Appendix [CJ The spectrum is symmetric un- 
der a change of sign in energy if we keep all energy levels 
into account. Andreev bound states arising from the gap 
edges move to the Fermi energy as Lp/ao increases^. 
There is level repulsion between Andreev bound state and 
oscillations of the energy levels with a period of order 10 
times £ f = 2tpao/Trh cx . These oscillations do not exist 
in Ref . [28[ Among all Andreev bound states obtained at 
a fixed Lp/ao some have a spectral weight much larger 
than the other. We keep only the levels having a spectral 
weight larger than a given cut-off. The evolution of the 
remaining Andreev bound states as a function of Lp/ao 
is then in agreement with Ref. |2^. 



In connection with the discussion in section llVB II we 
see that Andreev bound states near the Fermi energy 
occur only if the ferromagnetic chain is long enough, 
larger than approximately 10 times £s — 2isao/7rA. This 
should be contrasted with the LDOS in the infinite pla- 
nar limit discussed in section llVB II and with the model 
discussed in section IVII Bl 
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FIG. 11: Density of Andreev bound states (in arbitrary units) 
p(E/ A) as a function of reduced energy E/ A for Z = 1000 
and Z = 200 (a); Z = 40 and Z = 8 (b). We used the 
parameters A/e F = 0.02, A = 0.01, h eK /A = 0.5, d = A F /2, 
W ~ 15900Af where Xf = 2n/ltF is the Fermi wave length. 



B. Bogoliubov-de Gennes equations 

We consider now the model proposed by Cserti et ai— 
in which a 2D ferromagnetic dot with a rectangular shape 
of dimensions (d, W) is connected to a superconductor 
(see Fig.nj(b)). The superconductor has a width W but 
is infinite in the other direction [d! — +oo on Fig. ^| . Us- 
ing the solution of this modelSS based on Bogoliubov-de 
Gennes equations we generated the set {E n } of Andreev 
bound state energies. Like for the de Gennes-Saint James 
models the density of bound states is not equal to the 
LDOS in the superconductor—. The evolution of the den- 
sity of Andreev bound states as a function of the interface 
transparency is shown on Fig. 1111 The interface trans- 
parency is parametrized by the dimensionless coefficient 
Z, equal to the repulsive interface potential in units of 
the Fermi energjsi&. For large values of Z (correspond- 
ing to tunnel interfaces) we obtain an Andreev bound 
state miniband around E/Aq — /i cx /Ao- As the inter- 
face transparency decreases the miniband broadens and 
splits into two separate minibands (see Fig. 1111) . This is 
in a qualitative agreement with section lTV B II (see Fig.^J- 
(c) and Fig.01-(d)) where we obtain also Andreev bound 
state minibands that evolve inside the superconducting 
gap and can generate a large density of states at the 
Fermi energy for some values of the interface transparen- 
cies. 



VIII. CONCLUSIONS 



a minigap, which explains why the self-consistent super- 
conducting gap is weakened by the formation of an An- 
dreev bound state miniband near the Fermi energy. In- 
creasing temperature decreases the intensity of the An- 
dreev bound state peak in the LDOS which correlates 
with a reduction of the low temperature self-consistent 
superconducting gap. 

We found non monotonic temperature dependences 
of the superconducting gap with weak ferromagnets for 



Lp,Ls <C ^jr^jCs as well as for L$ 



d X) and 



Lf « &\ 



In the regime Ls <C £g and Lp > £^ 



we obtain non monotonic temperature dependences for 
some values of Lp/ do for strong ferromagnets. 

For the F/S/F trilayer in the parallel alignment we find 
pairs of Andreev bound states at opposite energies but 
there exists a minigap so that the self-consistent super- 
conducting gap decreases monotonically with tempera- 
ture. However at a fixed temperature the self-consistent 
superconducting gap is non monotonic as a function of 
the exchange field which is due to pairs of Andreev bound 
states at opposite energies in the regime /i ex /Ao ~ 1 <C 
t/A . 

Concerning the induced exchange field in the F/S/F 
trilayer in the parallel alignment we find that for weak 
ferromagnets the exchange field is in the same direction 
as in the ferromagnetic electrodes whereas it is in the op- 
posite direction for strong ferromagnets. The change of 
sign is visible both in the LDOS and in the magnetiza- 
tion profile. As the exchange field in the ferromagnets is 
increased we find a first order transition to the normal 
state, as for paramagnetically limited superconducting 
films in an applied magnetic field. 
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APPENDIX A: RESONANCES IN LATERAL 
CONFINEMENT 



In this Appendix we account for the differences be- 
tween the resonant and off-resonant F/S/F trilayers. The 
wave function in the lateral direction within a given elec- 
trode is described by the tight binding Hamiltonian 



We have presented a detailed analysis of F/S/F trilay- 
ers with weak ferromagnets based on microscopic Green's 
functions. The exchange field induced in the supercon- 
ductor in the regime h cx /Ao 3> 1 generates Zeeman split- 
ting of the LDOS the regime /i ox /A <C 1 is characterized 
by Andreev bound states. In the non perturbative regime 
t/ Aq ~ h cx / Aq ~ 1 the Andreev bound state miniband 
crosses the Fermi energy as the interface transparency 
is increased therefore competing with the formation of 



H = toY,{\ z + l )( z \ + \ z )( z + l \], (Al) 

z 

where to is the hopping between neighboring layers, and 
z/clq is an integer between to Lja^ — 1. The eigenstates 
with open boundary conditions are given by 
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and the energy is given by 
e n {L)=2t ± i 



nnao 



L 



«o 



(A3) 



with n = 1, L/a,Q. Let us consider a F/S/F trilayer 
with Lp/ag layers in the ferromagnets and Ls/ao layers 
in the superconductor. Quasiparticles in the supercon- 
ductor with a transverse quantum number ns can tunnel 
in the ferromagnets only if there exists an energy level 
with quantum number np close to resonance in the fer- 
romagnets, such that e ns (Ls) — e nF (Lp). If tunneling 
between the ferromagnets and superconductor is not res- 
onant then the F/S/F trilayer behaves as if the ferromag- 
nets were insulating. The lowest off-resonant values of 
(L s /a ,L F /a Q ) correspond to (L s / 'a , L F / a ) = (1,2), 
(1,4), (2,1), (2,3), (3,2), (3,4), (4,1), (4,2), (4,3). The 
lowest resonant values of (L5/00, Lp/ao) correspond to 
(Ls/ao, Lp/a ) = (1, 1), (1, 3), (2, 2), (3, 1), (3, 3). 



APPENDIX B: DISORDER IN THE F/S/F 
TRILAYER WITH ATOMIC THICKNESS 



where n a is the concentration of impurities and u a is the 
scattering potential. We suppose that u a is independent 
on wave vector. 

The fully dressed Green's function of the connected 
trilayer without vertex corrections is obtained through 



Get a — 



J - n( d ) f n [a >t - n [a > t un^.'h 
1 ya,a L a,ay a ,a L a,a ya,ct L a,byb,b ^b,a 



id). 



(B4) 

where G a , a stands for G QQ (k, ui), and gf) stands for 

-(d) 



g^f(k, (J), with i — a,a,b. 



2. Vertex corrections 

Lowest order vertex corrections correspond to pro- 
cesses in which a quasiparticle of the superconducting 
layer makes an excursion in one of the ferromagnetic lay- 
ers in between two scatterings on a given impurity in the 
superconducting layer. The fully dressed Green's func- 
tion is obtained through 



In this Appendix we give the Dyson equations of a 
F/S/F trilayer in the presence of disorder. 



1. Disorder in the superconducting and 
ferromagnetic layers 

We first neglect vertex corrections that are discussed 
in section ITY G 21 We replace the Green's functions g a , a 

of the superconducting layer by the Green's function gcf}a 
of the superconducting layer in the presence of disorder. 
To second order in disorder the Dyson equation takes the 
form 

5i d a( k > w ) = 9a,a(k,uj) + n a alg a , a (k, u) (Bl) 
dk' 



r« a (k-k')ff a , a (k',a;)« a (k'-k)M(k,a;), 



where n a is the concentration of impurities and the over- 
line denotes an averaging over disorder. Eq. I|B1|) is 
solved according to 



(d),i,i Ad), 1,2 

ya.a Ja, a 
Ad), 2,1 (d),2,2 

J a, a </q',q 



1 

V 



fl,2 
J a, a 
f,\ q 2,2 



(B2) 



yi,i 
^9 



where d = 



(/i;2) 2 and 2? = 1 + 2f a , a E f - 

(S/) 2 . 



M 2,2 

1,1 yl,l _ 2 : 2 y2,2 _ 
act, a g aa,a g 

The self-energies due to disorder take the form 



with £ d = E^E 2 ,' 2 



^9 

v 1.2 



E 2,2 



2 2 
1 a ll a (lQ 



(B3) 



G a , a = g { a % (B5) 
a {d) i a (d H G + a {d) i h a (d) L G 

ila,a l 'a,aya,a l 'a,aya,a i y a,ct v a,byb,b l 'b,a^'a,a 

dk' 



n a a Q g a 



(2n} 



-u a (d) t a (d) t a (d) u G 



23(d) / dk ' fi d '(d) f d Kd)r g/(«0 fi a 



where we used the same notation as for Eq. (|B4|) and 
9& = 9^(^M,u a = u a (k-k'). 



APPENDIX C: ALGORITHM FOR THE ID 
MODEL 

In this Appendix we detail the algorithm by which we 
calculate the energies of the Andreev bound states of the 
ID models. In the superconductor the spectral repre- 
sentation of the local propagator is obtained by summing 
Eq. JfJJl over all energy levels of the ID chain with open 
boundary conditions: 



1,1 1 \ 



L s /ao , 

2ao ST* • I x a + a o 
— y sin nil— 

Ls + ao 1 \ L s + a 

71—1 



(CI) 



x sin nn 



x/3 + ap 
Ls + a 



it: 



t'T, 



uj — E n — iij u> + E n — irj 



where a and (3 correspond to two sites in the ID chain at 
positions x a and xp. Similar expressions are obtained for 
the "22" and "12" components. In the ferromagnet the 
energy levels are given by = e n (Lp) — crh ex , where 
e n (L) is given by Eq. i|A3|l . The local Green's function 
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at the extremity of the ID ferromagnet is given by 

Mr \ 2a o • 2 ( n7Ta o \ 1 

LF + a °n~i \Lf + ooJ u)-e£' -in 

We denote by a a site in the superconductor, cho- 
sen far from the boundaries. At site a is connected 
the extremity "a" of the ferromagnetic chain. We note 
t = t a , a = t a , a - The fully dressed Green's function G]^ a 
of spin-up electrons at site a is deduced from the Dyson 
equations2&2£i2£. The spectral representation of G]^ a is 



obtained by evaluating numerically the position of the in- 
dividual energy levels w n and their spectral weights R n : 

LDOS is given by 

p&M - ~Im [<?&(<")] = £ ^'^(w ~ (C3) 

n 

The energy levels and spectral weights can be obtained 
without approximation to an arbitrary precision by using 
a dichotomy algorithm. 
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